Two-component approach for thermodynamic properties in diluted magnetic 

semiconductors 
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We examine the feasibility of a simple description of Mn ions in III-V diluted magnetic semiconduc- 
tors (DMSs) in terms of two species (components), motivated by the expectation that the Mn-hole 
exchange couplings are widely distributed, especially for low Mn concentrations. We find, using 
distributions indicated by recent numerical mean field studies, that the thermodynamic properties 
(magnetization, susceptibility, and specific heat) cannot be fit by a single coupling as in a homoge- 
neous model, but can be fit well by a two-component model with a temperature dependent number 
of "strongly" and "weakly" coupled spins. This suggests that a two-component description may be 
a minimal model for the interpretation of experimental measurements of thermodynamic quantities 
in III-V DMS systems. 



PACS: 75.50.Pp, 75.10.Lp, 75.40.-s, 75.30.Hx 
I. INTRODUCTION 

Diluted, magnetic semiconductors (DMSs) have been 
the focus of intense study recently due to their potential 
for use in novel devices making use of bath magnetic and 
conventional semiconductor propertiesJj'B The discovery 
of a magnetic transition temperature, Tc, of 110 K in a 
sample of Gai-ajMn^jAs for x — 0.053 has further spurred 
efforts to understand the origin and physical effects that 
influence the magnetic properties of these materials.E 

It is now reasonably well established that III-V systems 
such as Gai_a;Mn2;As with x = 0.01 — 0.07 are itinerant 
ferromagnets, in which the Mn ions play the dual roles of 
acceptor site and magnetic ion, and the itinerant carri- 
ers are holes, which-have an antiferromagnetic interaction 
with the Mn spinsuij The antiferromagnetic hole-Mn in- 
teraction leads to an effective ferromagnetic interaction 
between Mn spins, and gives rise to the ferromagnetic 
transition. One experimental fact that may be of impor- 
tance is that in these systems the number of holes is 
only a small fraction of the number of Mn dopants (or 
Mn spins), UMm implying that the system has low carrier 
density, and is heavily compensated. 

Theoretically, there have been several approaches to 
trying to calculate the thermodynamic and transport 
properties of these materials - one j^approach has been 
to look at the effect of spin wavesBllil whilst anothea; 
has been a mean field model including spin-orbit effects.u 
Both of these approaches leave out the disorder due to 
the random positions of the Mn ions in the sample. The 
effect of random positions has been considered in a nu- 
merical mean field theorylljllj for the low density phase 
of Gai_2;Mna;As, as well as in Monte Carlo simulations 
of the insulating phase of II- VI DMS, represepted by a 
Heisenberg model for the Mn and carrier spins .ll^'tj Both 
of these investigations show that the positional disorder 
gives rise to a distribution of exchange couplings between 
Mn ions and holes. Monte Carlo results have also been 
obtained for a model in which the hole-Mn coupling is 
assumed constant and leads to an effective Mn-Mn in- 



teraction where positional disorder is included.lll A re- 
cent Monte Carlfl-study of Gai_a;Mna;As using a kinetic- 
exchange modeljEZI has appeared whilst this work was be- 
ing written up. 

In this paper, we construct the simplest mean field 
model that attempts to capture the effects of disorder in 
the effective local fields at different Mn sites. This disor- 
der arises as a result of the different local potential for the 
carrier at different Mn sites. At one extreme we consider 
a simple model of compensation by antisite defects. In 
this model, each As antisite defect is viewed as capturing 
the holes of two neighboring Mn dopants, and providing 
an onsite potential that is significantly different from Mn 
sites that are far from such As antisites. This naturally 
leads to description of the Mn spins in terms of two dis- 
tinct species. This is clearly a caricature, since for posi- 
tionally random doping and As antisite defects, there will 
be a continuous distribution of onsite potentials rather 
than a bimodal one. However, as we show in the bulk of 
our paper, if the distribution is rather wide fei-j^ found 
in the mean-field study of Mn impurity bandal3'E3 ) , such 
a bimodal distribution provides a reasonable description 
of the thermodynamics, provided we allow the relative 
weights of the two species to be temperature dependent 
according to a simple rule, which occurs naturally in the 
analysis. 

We concentrate only on the carrier-spin exchange 
part of the. Hamiltonian, since a numerical mean-field 
treatmentllj shows that this term captures most of the 
condensation energy for the magnetic phase, and the car- 
rier kinetic energy changes only weakly with the onset of 
ferromagnetism. As stated earlier, in the regime of in- 
terest, the concentration of holes, n^, is considerably less 
than the concentration riMn of Mn ions. The fluctuations 
in the local carrier charge density (due to fluctuations in 
the occupation number of the impurity states around Mn 
sites with different surroundings) are then represented by 
a fluctuation of the effective exchange coupling between 
the Mn moment and the spin of the carrier. 

Our results can be summarized simply as follows : we 
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find that representing the hole-Mn antiferromagnetic ex- 
change coupUng by a single parameter is insufficient to 
capture the thermodynamic behaviour, such as the tem- 
perature dependence of the magnetization, susceptibihty 
and specific heat in the ferromagnetic phase, when dis- 
order is large. However, by using a model where there 
are two species of Mn ions, with different hole-Mn ex- 
change couplings, we are able to fit the magnetization 
and other thermodynamic parameters that we calculate 
from a distribution of couplings in a much more satis- 
factory manner. Furthermore, as argued below, such a 
two-component model captures the inadequacies of the 
homogeneous model in a manner that is qualitatively cor- 
rect for temperatures not too close to Tc- We emphasize 
that the label "component" used throughout this paper, is 
synonymous with "species " or "type ", and has no relation 
to the number of components of the spin itself, which in 
this work refers to a vector in three dimensions. 

The idea of using a two-component model is similar in 
spirit to models developed to understand the magnetic 
behaviour of doped non-magnetic semiconductors such 
as phosphorus doped silicon, motivated by the observa- 
tion that exchange couplings in such systerns .should be 
distributed over many orders of mag nitudclalia In Mn 
doped DMS systems, disorder appears to lead to a simi- 
lar situation of widely distributed couplings .t3 Whilst not 
explicitly included in Refs. |l^,|l^ the large compensation, 
apparently from As antisite defects, would enhance disor- 
der, and hence broaden the distribution of effective cou- 
plings. Therefore a two component approach appears to 
be a natural approximation that is qualitatively correct, 
and may be quantitatively adequate for many purposes. 



II. MEAN FIELD MODEL 



A. Hamiltonian 

The DMS system we study consists of magnetic ions 
(Mn) coupled to charge carriers. In the case of II- VI 
semiconductors (e.g. ZnSe), the carrier is provided by a 
second dopant (such as P); however, in III-V semiconduc- 
tors (e.g. GaAs), Mn, being a divalent atom substituting 
on a trivalent (Ga) site, provides a hole in addition to the 
spin. In this paper, we use the example of Gai_a:Mnj;As 
for concreteness. 

For Mn concentrations of interest (a fraction of a per- 
cent to a few percent), and hole concentrations of about 
10%-20% per Mn impuritv- [the system is near a metal- 
insulator transition (MIT)J3'E3 This imphcs that the hole 
wavefunctions are filamentary, with a multifractal struc- 
ture on length scales that determine the magnetic be- 
haviour of the system, which is quite distinct from the 
homogeneous structure for plane wave (or Bloch wave) 
states characteristic of periodic systems. Thus, each hole 
interacts with many Mn spins, depending on the ampli- 
tude of its wavefunction at various Mn sites, as well as the 



envelope function characterizing the hole. Since each hole 
interacts with a large number of spins the net exchange 
fields felt by different holes are similar (i.e., the fluctua- 
tions are not that large). However, the hole concentra- 
tion is considerably (a factor of 5-10) smaller than the 
Mn concentration; therefore, each Mn spin experiences a 
rather different exchange field due to the few holes that 
have significant amplitude at that site (or nearby sites, 
via the tail of the envelope function). Hence the fluc- 
tuations in the local exchange field at different Mn sites 
cannot be ignored, since the fields are being produced 
predominantly by just a few holes. This-asymmetry has 
been documented in a numerical studyllj and forms the 
basis of our simple phenomenological scheme. These fluc- 
tuations are of paramount importance in the insulating 
phase and at the MIT. If experience with conventional 
doped semiconductors is any guid&™ .-tiiey are likely to 
persist well into the metallic phasala^t^ and cannot rea- 
sonably be ignored in any theory of DMS ferromagnetism 
up to dopant densities well above (a factor of 3-5) the 
dopant density for the MIT. The above discussion raises 
questions about the applicability of the results of studies 
which are based on homogeneous mean-field models, as 
well as those based on perturbative treatments of the Mn 
spin system (such as RKKY exchange), given the large 
ratio of spin to carrier density. 

The dilutcncss of the carrier system leads us to neglect 
hole- hole interactions. Direct Mn-Mn interactions are 
also ignored because they are extremely short range due 
to the atomic nature of the Mn 3d-orbitals responsible 
for the Mp_spin. Furthermore, the numerical mean field 
treatmentli3 shows that the energy gain due to the onset 
of ferromagnetic ordering coming from the exchange term 
in the Hamiltonian is much larger than the change in 
carrier kinetic energy. Ignoring the kinetic energy is not 
expected to lead to any major qualitative differences. In 
fact, in our results it does not lead to even quantitative 
differences in the magnetization for temperatures below 
about 0.6 Tc (as shown in Section HI), while there are 
quantitative differences nearer the ordering temperature, 
as the exchange energy and kinetic energy variations with 
T become comparable in magnitude. Consequently, in 
keeping with the philosophy of finding a minimal model 
description, we consider only the exchange term of the 
Hamiltonian in this paper. This Hamiltonian takes a 
form similar to that studied by Konig et aln with no 
kinetic term: 



J d^r^ J,„(r)s„ • S,(r) 



5/iBB.^S,(r)+g>BB.^s„(r) , (1) 



where a labels the holes, i labels the Mn spins, is the 
spin for the a*'' hole and Si (r) is the Mn spin centered on 
site i. The g-factors of the hole and Mn spins are labelled 
by g* and g respectively, /is is the Bohr magneton and 
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B is an external magnetic field, which we shall assume 
to be zero unless otherwise stated. The overlap integral 
for the a*'* hole with the spin centered on site i is Jia- 
We assume that the d-electrons that give rise to the Mn 
spins are localized in comparison to the holes and treat 
the Mn spins as having delta function spatial dependence, 
i.e. Si(r) — SiS'^^^r — R^), which, after integrating out 
the delta functions, leads to the Hamiltonian 



H — ^ ^ Jia{'R-i)Sa ■ Si 



(2) 



where JiaCRi) — Jo|0Q(Ri)|'^, with (f>a{'r) the wavefunc- 
tion of the a*'* hole and Jq the microscopic exchange con- 
stant. We can write J2a '^ia{'^i)^a = h(Ri), and treat 
the effect of the magnetization of the hole spins as creat- 
ing an effective field at each Mn site, so our Hamiltonian 
is 



7^-^h(R,)-S, 



(3) 



In a recent mean-field studj^l of a tight-binding model of 
the impurity band arising from holes on Mn sites coupled 
to Mn spins, the randomness of the Mn sites is explicitly 
taken into account. That work involves a numerical im- 
plementation of the self-consistent mean field equations, 
in which the effective fields are self-consistently calcu- 
lated at each Mn site for each temperature, T. (The 
mean-field in such treatments refers only to temporal av- 
eraging over the local environment, and not positional 
averaging). A ferromagnetic phase is found below a crit- 
ical temparture, with a distribution of effective fields, 
P{h), which is temperature dependent. The numerical 
mean field model explicitly includes the itinerant nature 
of the holes. We use the distribution of effective fields 
from the self-consistent mean-field calculations. Hence, 
even though the models constructed here are in the form 
of pure exchange, the wavefunctions used to calculate 
effective couplings are those of itinerant electrons. Con- 
sequently, that physics is implicit in these models. In 
a mean field description of carrier moments, the field is 
simply a product of the mean moment (which we take 
to be along the z direction, of magnitude s^) and an 
effective exchange interaction J which varies from site 
to site. Thus, we can obtain from the numerical study, 
an effective distribution of exchange couplings -P(J), and 
our Hamiltonian becomes an integral over the distribu- 
tion P{J): 



istributii 



dJP{J)JS^js' 



(4) 



Details of the derivation are given in Appendix 

In a fully self-consistent scheme (such as Ref. |12| ) the 
distribution P( J) depends on the temperature T due to 
the temperature dependence of the hole amplitude dis- 
tribution. However, as we show below, the distribution 



Po{J) calculated at T = works very well for tempera- 
tures up to T = 0.6 Tc- Hence, in keeping with our search 
for a simple model, we will use Po{J), and consider the 
T-dependence of P{J) only where necessary. 



1. Spin and carrier magnetization 

We find a self-consistent solution for the average Mn 
spin of 



(5}) = ~SBs{f3aj), 
where S' = 5/2 is the spin of each Mn ion, 

1 



" 25 i — 2 — ^ 



■ coth , 
25 V2 



(5) 



(6) 



is the Brillouin function for spins with magnitude S and 
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We also need to take into account the fact that the holes 
experience an effective field from the Mn spins, which 
leads to the following self-consistency condition 



where 



a 



P JO 



dJ P{J)J (S^j) , 



(8) 



(9) 



p is the ratio nh/nyin and j is the effective hole spin 
(j — 3/2 in real syst|e*n*j|. To compare our results with 
the numerical studyjij't3 we use the Brillouin function 
with j — 1/2 for the hole spin and replace Jo by 3 Jo 
so that in effect the hole spins take values of ±3/2. We 
emphasize however, that we have found that treating the 
spin as a classical object (i.e. using a Langevin hmction) 
or using Brillouin functions with j — 1/2 or j = 3/2 
does not affect the basic picture we have here, or the 
characteristic shape of the curves, other than in raising 
or lowering Tc. The value of Tc is given by 



(10) 



where denotes an average with respect to P{J) and 
(. . .) denotes a thermodynamic average. Note that whilst 
p enters Equation ( p^ ) explicitly, there is also an implicit 
dependence, in that p affects the distribution of hole-Mn 
exchange couplings P{J) and hence J^. 



2. Susceptibility and Specific Heat 

We next calculate the magnetic susceptibility and spe- 
cific heat at zero field. The susceptibility per unit volume 
is 
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X = X + 



dJ Xj, 



(11) 



where x* is the contribution to the susceptibihty due to 
the holes and xj is the contribution to the susceptibihty 
from Mn spins coupled to holes with exchange J. The 
individual expressions for the susceptibilities are 



X 



XJ = gnBP{J)nMnf3Gj ^g^B 
where 



dJP{J).PGj 

1 J * 
-X 



nh g*tJ.B 



(12) 



(13) 



G* 



(14) 
(15) 



We use the expressions for the susceptibility given in 
equations ( |l2| ) and ( |l3| ) with g — 2, and assmne g* = 2 
for our numerical calculations in section H.Eil 



To obtain the specific heat we have to be careful since 
we calculate with a Hamiltonian that has temperature 
dependent coefficients and the concentrations of strongly 
and weakly coupled spins may also depend on temper- 
ature. We use the mean field free energy to obtain the 
specific heat through the relation 



dU 
'dT' 



(16) 



where we note that the energy per Mn spin is 



1. Single component model 

In a single component model of Mn spins, where each 
spin is coupled in the same way to the carrier spins, the 
form for P{ J) is just a delta function: 



P{J) ^ S{J - Ji), 



(18) 



where Ji is the exchange coupling. In that case, the 
formulae for thermodynamic quantities derived for a dis- 
tribution in the previous section lead to 



{s^} =-jB,{(3a*), 



(19) 
(20) 



with a* = {l/p)J (S^) and ai = Ji (s^) . The suscepti- 
bility per unit volume is 



X = X +Xi, 



and 



{g*^lB)'nhPG* (l-^|jiGi) 
^ " 1 - ^(3^G*JlGi ' 

Xi = (5MB)"-Mn/3Gi (guB - — -r— X" 
V nh g*HB 



(21) 

(22) 
(23) 



where the G functions have the same meaning as previ- 
ously. Finally we calculate the specific heat again with 
the derivative of energy with respect to temperature, but 
in this case the energy per Mn spin is 



U 



(24) 



U 



Ni 



Mn 



dJ P{J)J (S^j) (s^) , 



(17) 



in our mean field theory, where A^Mn is the total number 
of Mn spins. 



B. Single and two component models 

For a system without disorder, the exchange between 
carriers and spins can be characterized by a single cou- 
pling, as has been done, e.g., in Ref. ^ To demonstrate 
how these models are inadequate for the system under 
consideration, we consider two simplified models to sub- 
stitute for the hierarchy of couplings implied by the dis- 
tribution Po(J). The first has a single coupling parame- 
ter, while the second incorporates the idea of "strongly" 
and "weakly" coupled spins in terms of two coupling pa- 
rameters. In both cases, the parameters are determined 
from the distribution Po(>^)- We compare the results of 
these two models with that obtained from the distribu- 
tion Po{J) in section |l| . 



2. Two component model 

As we show in the next section, the calculated distri- 
bution of exchange couplings are quite broad, and cover 
many orders of magnitude for parameter values of inter- 
est. This motivates us to study a model that is the next 
simplest after the single coupling model, namely the two 
component model, with an exchange distribution: 

P(J) = — (5(J- Ji)-K— (5(J- J2). (25) 



The physical motivation of the above distribution is to 
divide the Mn spins into two types, one which is strongly 
coupled to the carrier spins (Ji), with an effective con- 
centration ni, and the other which is weakly coupled to 
the carrier spins (J2) with a concentration 712. Since the 
only energy scale characterizing the thermodynamics is 
the temperature T, we would expect it to play an impor- 
tant role in determining both the coupling constants Ji 
and J2, and the concentrations of strongly and weakly 
coupled spins. Consequently, we expect that the best fit 
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to the curves for a distribution of exchanges will be ob- 
tained when Ji and rii are temperature dependent. This 
temperature dependence of parameters should be viewed 
in the same spirit as in a variational fit to free ener- 
gies of actualj-lT-independent) Hamiltonians by model 
Hamiltonians.ES 

We obtain a self consistent mean field solution to this 
model of 



{s^) =-jB,{pa*), 

where 

with a = 1 or 2, and 

rih rih 



(26) 
(27) 

(28) 
(29) 



Using the notation introduced above, the susceptibility 
per unit volume is 



X = X +Xi +X2, 



(30) 



where x* is the contribution to the susceptibility due to 
the holes and xi,2 are the contributions from the two 
species respectively - the expressions for the susceptibil- 
ities are 



1 - 4^/3JiGi - A^P,hG2 



lf32Q*rjl^j2Q ^ ^J2G2) 



Xa = gi^BTLaPGa ( 5MS - — ^^X' 

nh g*lJ.B 



(31) 
(32) 



where the G functions are as defined in Equations (|lj) 
and ( |l5|) . The energy per Mn spin (from which we deter- 
mine the specific heat) is 

JL ^ J^j^ (,^) ^st) + {s^} (S^) ■ (33) 
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FIG. 1. The distribution Po(logj^Q J) calculated for p = 0.1 
and X — 0.01 in the numerical mean field model of Ref. 13 at 
T = 0. 



III. NUMERICAL TESTS OF THE ONE AND 
TWO COMPONENT MODELS 

We now compare results obtained using our simple ap- 
proximations with the full numerical mean field calcula- 
tions. We concentrate on two cases, where the Mn con- 
centration UMn leads to a fractional occupancy of the Ga 
sublattice by Mn ions of a: = 0.01 and 0.02. For both 
cases, we consider a ratio of hole concentration rih equal 
to 10% of the Mn concentration, i.e. n/i/jiMn "r-P — 0.1. 
The distribution Pq{J) calculated numericallyll3 for the 
X — 0.01 case at T = 0, using a Bohr radius for the hole 

equal to 7.8 A,@ is shown in Fig. 0. As can be seen, 
it spans almost three orders of magnitude, and for this 
density, consists of two peaks. The higher peak is due 
to sites where the exchange interaction is dominated by 
a single hole that has a high probability of being on the 
Mn site in question, whilst the lower peak is found to be 
due to sites that have practically no amplitude for a hole, 
but whose exchange field is coming from holes on nearby 
sites. As X and p are changed, the relative weights in 
the two peaks changes, as does the total width of the 
distribution. However, the inferences made for this con- 
centration remain valid for higher values of x (we have 
checked explicitly the cases x = 0.02 and x = 0.03). 



1 ' 






■i \ Pt,(J) 
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Numerical mean field 


-..Pom 



T/Jo 

FIG. 2. Illustration of how the choice of P(J) at different 
temperatures for x — 0.01 influences the fit to the numerical 
mean field magnetization curve. 

Figure |2| plots the average Mn spin, Suni T), c alculated 
using the self-consistent solution of section [I A with the 
distribution Po{J) (Fig. |l|) (dashed line) and the full 
numerical mean field result (solid line) against tempera- 
ture. As can be seen clearly, the curve using the T = 
distribution works very well (errors less than 1%) until 
about 60 % of Tc for the numerical mean field model. To 
fit the numerical results properly for higher T, one must 
allow for the distribution of local fields to deviate from 
the T = distribution as the polarization of the holes be- 
gins to fluctuate strongly as the transition is approached. 
In figure ^, we indicate with crosses the SMn{T) curve 
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obtained using the distribution Pt^{J) (determined for 
T = Tj,)- Clearly, an interpolation of distributions at the 
two extremes {T — Q and T — Tel— will be adequate to 
reproduce the full numerical curve.ES 

Our goal is, however, to simplify the description of 
the distribution of exchange fields in terms of a few cou- 
plings. To this end, we consider the case with the fixed 
T-independent P{J) shown in Fig. |l|, and attempt to fit 
the thermodynamic properties for that case using a one- 
or two- component model of Mn spins. To fit the full 
numerical results for temperatures up to (for the nu- 
merical mean field model) , we would use a similar scheme 
with the appropriate PJJ) which best fits those numeri- 
cal results. In Figures |3| to ^, Tc — 0.84 Jq refers to that 
determined from Po{J) as shown in Figure |2[ 



ni 




T/T, 



FIG. 3. Comparison of Mn magnetization as a function 
of temperature for the numerically derived distribution P{J) 
(solid line), the single coupling J = J (labelled, crosses) model 
and the two component model with linear averaging (dashed 
line) where 7 — 0.6 and x = 0.01. 

For the single coupling parameter model, we replace 
the distribution P{J) by the mean coupling 



.h =J = 



dJ JP{ J), 



(34) 



which foe, X = 0.01 when we use Po{J) gives J — 
0.161 JqE3 Figure ^ shows the average Mn spin SuniT) 
(dotted curves) calculated using this model for x = 0.01. 
The results are clearly seen to be a poor fit to the results 
obtained for the full distrii).ution (solid lines), over most 
of the temperature range.Es 

In the two-component model, the parameters are de- 
termined using the following scheme. First, since the 
temperature is the only thermal energy scale in the prob- 
lem that can be used to define "strongly" and "weakly" 
coupled spins, we define a cut-off coupling = jT . All 
spins that have couplings below Jc are weakly coupled, 
whilst those with couplings greater than Jc are strongly 
coupled. Thus the concentrations ni and 712 of the two 
sets of spins are given by the relative fractions of spins 
with couplings above and below Jc, i.e.: 



"•Mn J J, 

n2 _ ^ _ ni 



dJP{J), 



"•Mn "Mn 



(35) 
(36) 



On the other hand, the couplings Ji and J2 are taken to 
be the averages over the two populations 



Ji = J dJJP{J), 
J2 = J dJJP{J). 



(37) 
(38) 



In this scheme, the only adjustable parameter is the di- 
mensionless parameter 7, which determines Jc and is cho- 
sen to give the best overall fit to SuniT). Using Pq{J) 
yields 7 = 0.6. Since most thermodynamic functions 
depend exponentially on the ratio J/T, we expect this 
scheme to work especially well when the distributions 
are broad, and the demarcation betweea "strongly" and 
"weakly" coupled spins becomes sharp £3 For narrow dis- 
tributions, on the other hand, this scheme essentially re- 
duces to a single coupling, which should be adequate for 
most purposes. 

The curve of Mn spin versus temperature SMniT), ob- 
tained using the two-component model, is shown in Fig- 
ure ^ as a dashed line. In contrast to the single coupling 
model the curve has the same shape as for the full dis- 
tribution and provides a much better quantitative fit. It 
might appear that using two couplings to approximate the 
distribution shown in Fig. [J works well because of the 
double peaked nature of the distribution. We wish to em- 
phasize that this is not the case; the main reason actually 
appears to be the large width of the distribution Po{J). In 
particular, if the upper peak in Po{J) is removed, then the 
magnetization curve from the modified P{J) still needs a 
two- component model to provide an adequate fit and a 
one- coupling parameter fit works barely better than for 
the distribution shown in Fig. [J. 

We now turn to other thermodynamic quantities, the 
susceptibility and specific heat; these quantities are 
shown for x — 0.01 and were calculated for a distribu- 
tion, as well as the single and two component models, in 
Section ||. Figure |^ plots the susceptibility for the mean 
field model with the distribution P{J), the two compo- 
nent model with 7 = 0.6 and the single component model 
for X — 0.01. Because the susceptibility is a higher order 
derivative, the two component model is not quite as accu- 
rate as for the magnetization. Nevertheless, there is good 
agreement with the results for the full distribution on a 
semiquantitative level down to 10 % of Tc as obtained 
from Po{J), whilst the single coupling model bears little 
resemblance to the distribution. 

Similarly, the curves for specific heat (Fig. ||) as a func- 
tion of temperature show good agreement between the 
two-component model and the distribution (in fact, bet- 
ter than for the susceptibility) for temperatures greater 
than 0.1 Tc. In contrast, there is a strong quantitative 
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discrepancy between the single coupling model and the 
distribution. The Schottky type anomaly is broadened 
out considerably for the distribution as well as the two- 
component model; sirnilfir broadening is present in the 
full numerical solutionJla 




0.01 0.1 1 



T/Te 

FIG. 4. Comparison of the susceptibility as a function of 
temperature for x = 0.01 calculated for the distribution P{J) 
(solid line) , the single coupling J = J model (crosses) and the 
two component model (dashed line) with 7 = 0.6. Note that 
Tc is that for Po{J). 

Figure ^ apparently shows a divergence in the suscepti- 
bility at low temperatures - this is not a divergence, but 
a peak at very low temperatures, due to the considerable 
population of Mn spins with very small local fields. 

1 I . . — . — ■ ■ — 
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FIG. 5. Comparison of the specific heat as a function of 
temperature for x = 0.01 calculated for the distribution P{J) 
(solid line) , the single coupling J = J model (crosses) and the 
two component model (dashed line) with 7 = 0.6. Note that 
Tc is that for Po{J). 

Finally, to understand better the behaviour of the two- 
coupling model, we show the temperature dependence of 
the couplings Ji, J2 and Jc and also the temperature 
dependence of the ratio ni/nMn in Figures I and re- 
spectively. 
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FIG. 6. The couplings Ji and J2 for x = 0.01, shown with 
Jc ~ 0.6 T as functions of temperature where Ji and J2 are 
obtained in the linear averaging scheme. Note that Tc is that 
for Po{J). 

As expected from the physically motivated criterion 
described earlier in this section, at low temperatures 
nearly all of the spins are strongly coupled; however, 
as the temperature rises the number of strongly coupled 
spins decreases sharply, as the strength of the coupling 
required for a spin to be strongly coupled increases with 
T. It is this essential aspect of a broad distribution of 
couplings that allows the two-component model to prop- 
erly capture this behaviour. 
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FIG. 7. The ratio ni/nmn as a function of temperature for 
X = 0.01 using the model in which Jc — 0.6 T where Ji and 
J2 are obtained in the linear averaging scheme. Note that Tc 
is that for Po{J). 



IV. DISCUSSION AND CONCLUSIONS 

Whilst other authors have acknowledged that dis- 
order plajiSi-soHie role in the physics pf.|4JI-V DMS 
sy stems, □'tiHl30 the resiJ±s-of mean-fieldliSH as well as 
Monte Carlo simulationsEj'EJ on models with a random 
distribution of dopants near the metal-insulator transi- 
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tion suggest that there should be strong effects on the 
thermodynamic properties of DMS systems, especially at 
low temperatures. The tendency in experimental fits, on 
the other hand, has been to neglect disorder effects en- 
tirely and produce fits characteristic of homogeneous sys- 
tems. Thus, magnetization curves obtained from jtijans- 
port measurements in a sample with x = 0.035,c3 are 
fitted with Brillouin functions to extract a single ex- 
change coupling. It is important to ask whether the 
effects described here have been seen or have the po- 
tential to be seen experimentally. Significant deviations 
from a Brillouia. ci.irye appear to be seen in a number 
of experimentsj3'EIrc3 in contrast to Ref. The Bril- 
louin function-like behaviour appears to arise in trans- 
port, whilst the deviations are seen in SQUID measure- 
ments of the magnetization. One explanation for this 
may be that the transport measurement mainly sam- 
ples Mn spins that are strongly coupled to holes, which 
can give a Brillouin function as described below. How- 
ever, the SQUID measurement samples all Mn moments 
equally, and the magnetization inferred is from both 
weakly and strongly coupled Mn spins (as is calculated 
here). Interestingly, a measurement of the magnetiza- 
tion against magnetic field at a temperature of 2 K (Tc 
was 37 K) in a sample with x = 0.02 showed that the 
magnetization was roughly only 25 % of its saturation 
value.E3 This is consistent with our picture here that due 
to a broad distribution of exchange couplings there are 
significant numbers of spins that are not polarized, even 
at low temperatures. In studies of insulating samples it 
has been observed that the saturation moment was con- 
sisten^wil,h only about 40% to 50% participation of Mn 
spinsEJ^El 



2.5 
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FIG. 8. Illustration of how a single Brillouin function can 
approximate part of the magnetization curve of a distribution 
{x — 0.01). The parameters chosen for the single Brillouin 
function were that 28 % of the Mn spins were coupled to 
holes with J = 0.57 Jo, whilst the rest were uncoupled. 

Figure ^ illustrates how it is possible to have the mag- 
netization from the distribution Po{J) fitted by a single 
Brillouin function over some temperature range. The pa- 



rameters used to obtain the curve shown were to assume 
that only 28 % of the Mn spins contribute to the magne- 
tization and that the coupling is J = 0.57 Jq. Also, in the 
large x insulating phase, there has been an observation 
of a magnetization curve that could not be fitted, with a 
Brillouin function assuming only a single coupling .lIjH In 
that work this was ascribed to multiple exchange mech- 
anisms rather than disorder in the local fields as we sug- 
gest here. We believe similar effects will persist on the 
metallic side of the metal-insulator transition, and the 
assumption of a single exchange coupling used in experi- 
mental fits will be valid only deep in the metallic phase. 




B = oT ; I 

B^5T o 



T/K 

FIG. 9. Two component fit to SQUID magnetization curve 
at zero field and B = 5 T for a sample with x = 0.053 and 
p — 0.3 from Ref. 27. The experimental data (• for B = T, 
and o for B = 5 T) is fitted very well with a two component 
model with Ji — 47.5 K, J2 = 7.5 K and ni = 0.41nA/n (solid 
bold fine for B = T, dashed bold line for B = 5 T). The 
light curves are for a single coupling and J = 31 K (solid B 
= T and dashed B = 5 T). 

To illustrate this point, we show data for the rema- 
nent magnetization and field dependent magnetization 
from a SQUID measurement (Ref. |2^) for a sample with 
X = 0.053 and p = 0.30 that has been fit with j — 3/2 
at zero field and at i? = 5 T. (Note that the only change 
required to include a magnetic field in the formalism in 
Section H, is a — > a — gfisB). Since we obviously do not 
know P{J) in this case, and we expect it to be less broad 
in the metallic phase, we have taken temperature inde- 
pendent Ji , J2 and ni , and compared the fit to one with 
a single J that reproduces Tc correctly. Using Ji = 47.5 
K, J2 = 7.5 K and ni — 0.41njv/n we obtain the fits 
shown in Fig. |^ as bold solid and dashed lines, which 
are in close accord with the data at both magnetic field 
strengths, whereas a fit with a single J = 31 K (light and 
solid dashed lines) is clearly inadequate. Whilst a single 
fit of this form does not necessarily imply our model, 
it would clearly be of use to determine if a such a two 
component behaviour is universally seen, by fitting ex- 
perimental data for M{H, T) at several fields and tem- 
peratures. 



8 



In conclusion, we have found that we can approxi- 
mate the results of the numerical mean field treatment 
of a kinetic- exchange model of Mn dopants in GaAs, 
which properly treats the positional disorder in the alloy 
system, by a simple exchange model with temperature- 
dependent effective couplings between carriers and two 
different species of Mn spins. Such a model emerges nat- 
urally because of the distribution of couplings that are 
a consequence of the positional disorder of the Mn ions. 
Since the ultimate inputs in the simple model are the 
couplings Ji and concentrations Tij of each species, by 
parametrizing the behavior shown in Figures 7 and 8 in 
terms of a single parameter (or perhaps two), experimen- 
tal data can be used to yield information about the true 
nature of coupling distributions in actual DMS materi- 
als. We hope that with more detailed measurements, 
both of bulk magnetic and thermodynamic quantities, as 
well as those from local probes, it may be possible to 
parametrize the couplings accurately enough to provide 
a quantitative fit to the other thermodynamic properties 
at various concentrations in all regions of the phase dia- 
gram, in the quantitative manner possible for crystalline 
systems with translational symmetry. 
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APPENDIX A: DERIVATION OF INTEGRAL 
OVER EFFECTIVE COUPLINGS 

The Hamiltonian (Eq. (2)) may be rewritten as 

H = Jo^|0a(R.)|'s„-S„ (Al) 

where we have used the definition Jict(Ri) = Jo|0c<(Ri)p. 
We assume that the amplitude of the a**^ hole at site i 
may be written as the sum over Mn sites of the product 
of the amplitude for a hole to be at a given site times the 
amplitude from the local atomic orbital ipir) (note that 
ip{r) oc e^*"/"^). Hence 

\M^r)? = - R,)l', (A2) 



where we are in effect ignoring any quantum interference 
terms. This allows us to rewrite the Hamiltonian as 

n = Jo^p«(Rj)l^(R« - RjOl'sa • S.. (A3) 



One of our mean field assumptions (based cp. the idea 
that each hole interacts with many Mn spinsES), is that 
each hole behaves identically - in which case, (sq) = 
(s) independent of a. This also means that Pq(Rj) = 
(l/iV/i)p(Rj), where p(Rj) is the amplitude for finding 
any hole at site i. 
If we define 



J, = Jo^p(Rj)IV'(R» -Rj)!', 



(A4) 



then the effective exchange field at site i is Ji and the 
Hamiltonian takes the form 



To convert this sum into an integral, define 



P 



(A5) 



(A6) 



In mean field, (Sj) depends only on Jj, hence we can 
relabel Si — Sj and the Hamiltonian takes the form of 
Eq. (4) when we use the identity 



dJP(J) = 1. 



(A7) 
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